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Abstract Superfluid "^He at a finite temperature is described by the two-fluid 
model with the normal fluid component and the superfluid component. We formu- 
late the two-fluid model for two-component BECs, namely two-superfluid model, 
starting from the coupled Gross-Pitaevskii equations. The two-superfluid model 
well corresponds to the two-fluid model in superfluid ^He. In a special condition, 
the two sound modes in the two-superfluid model behave like first and second 
sounds in the two-fluid model of superfluid ^He. 
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1 Introduction 

Superfluid "*He has been thoroughly studied theoretically and experimentally in 
the field of low temperature physics since Kapitsa discovered superfluidity of ^He 
below the transition temperature ' . Tisza^ and Landau^ succeeded in understand- 
ing the superfluidity of "*He with introducing the two-fluid model, which states 
that the system consists of normal fluid and superfluid being independent of each 
other and is described by 



P„(^ + K-VK) =-^VP-p„aVT + r)ny^v„, (1) 

p,(^ + {v, ■ V)t;,) = -^VP + p,aVr. (2) 

Here p„ and d„ are density and velocity of normal fluid, and Ps and are those 
of superfluid. a is entropy per unit mass of the normal fluid, p = p„ + is total 
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density and t]„ is the coefficient of viscosity of the normal fluid. The pressure gra- 
dient VP runs both components in the same direction and the thermal gradient VT 
does in the opposite direction. Thermal counterflow driven by a thermal gradient 
is characteristic of superfluid "^He. When the relative velocity between two com- 
ponents is large, they become dependent through the mutual friction Fs„, which 
is added to Eqs. (1) and (2)'*. Other formulations for superfluid ^He are derived 
from the conservation law and equations of motion of mass density and entropy 
density. The hydrodynamic equations^ are 



dp 
dt 



= -Vj, (3) 



^ = -^V.K-..), (4) 



di 
dt 



VP, (5) 



l-{v„-Vs) = -^VT, (6) 

dt Pn 

where j = PnV„ + psVs- These equations yield the wave equation of mass density 
and entropy density, which leads to first sound and second sound. First sound 
is a mode of oscillation of total density and it exists generally in a usual fluid. 
While, second sound is a characteristic mode in superfluid ^He, in which entropy 
oscillates without oscillating total density, not existing in a usual fluid. 

An atomic BEC is one of the most important subjects in modern physics. Es- 
pecially, two-component BECs are known to create various exotic structure of 
quantized vortices^ and cause some characteristic hydrodynamic instability such 
as Kelvin-Helmholtz instability^ and Rayleigh-Taylor instability^. In another pa- 
per, we investigate counterflow in two-component BECs which has many analo- 
gies with thermal counterflow in superfluid ^He. For example, when the relative 
velocity exceeds a critical value, the counterflow becomes unstable and quantum 
turbulence appears like in thermal counterflow. In this work, we describe two- 
component BECs following the two-fluid model of superfluid "^He and obtain four 
equations similar to Eqs. (3)-(6). We derive two sound modes from these ele- 
mentary equations. Two sound modes in two-component BECs were obtained by 
some other works ^, but we have them correspond to first and second sounds. This 
is the main point of this work. Thus we can expect to improve interactive studies 
in superfluid "^He and two-component BECs with investigating common features 
between these. 



2 The two-fluid model in two-component BECs 

We consider binary mixture of BECs described by the wave functions Wj = ^Jnje"^j 
in the mean-field approximation at T = K, where the index j refers to each com- 
ponent j {j = 1,2). The wave functions are governed by the coupled Gross- 
Pitaevskii (GP) equations^S, 
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ih^^ = -;^V2'ft + y(r)f2 + g22|f2pf2 + gl2|fiP'ft, (8) 

where mj is particle mass associated with the species, gjj is intracomponent in- 
teraction and gi2 is intercomponent interaction. We insert Wj into Eqs. (7) and (8) 
and obtain the hydrodynamic equations, 

^ = -V(pi^i), (9) 

dpi 
dt 



-v-(p2t;2), (10) 
-Piv( 



^ mf mim2 

-P2V( 



2yP2 ^ '"2 

^fp2-^Pi). (12) 
2 mim2 / 

where pj — nijUj is mass density and Vj = ^^^j is superfluid velocity. Equations 

(9) and (10) are equations of continuity for pj and Eqs. (11) and (12) are quasi- 
Euler equations for the superfluid velocity. 

We will derive equations similar to Eqs. (1) and (2) to reveal correspondence 
between superfluid "^He and two-component BECs. Here we consider a uniform 
system and apply the long-wavelength approximation, so the potential term and 
the quantum pressure term in Eqs. (11) and (12) are neglected. The pressure of the 
whole system is 

p_ gllPl I g22p| ^ gl2PlP2 
2m J 2m2 niimi ' 

since P — —dE/dV. Then Eqs. (11) and (12) turn into 

Pi(^ + (i;i-V)t;i)=-ivP-ivf, (13) 

p2{^ + iv2-V}V2)=-^VP+^Vf, (14) 

with 

Vf = |i^Vp2 - I^Vpf + ^l^(piVp2 - P2VP1). 
2OT[ 2m2 mim2 

It is impossible to describe the right hand side by gradient of some scalar potential 
because pi and P2 are spatially dependent, but we represent it by Vf in order to 
emphasize the correspondence to Vr in Eqs. (1) and (2). From Eqs. (13) and (14), 
we can find that two-component BECs are driven by two terms. The pressure 
gradient VP runs both components in the same direction, while Vf runs them 
oppositely. This nature is just the same as one of the two-fluid model in superfltiid 
^He. 
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3 First sound and second sound in two-component BECs 



In this section, we will derive two sound modes from the four elementary equa- 
tions in two-component BECs and let them correspond to first and second sounds. 
Here we assume that superfluid velocities vj are small and the non-linear terms 
are neglected. By making Eq. (9) ± Eq. (10) and Eq. (1 1) ± Eq. (12) we obtain 

^p+ = -V.j+, (15) 
%_ = -V-j_, (16) 



dt 

d_ 



where p± = pi ± p2 and j± = piv\ ± p2V2- Because first sound means oscillation 
of p with two components in phase and second sound means oscillation of a with 
them out of phase, we expect that oscillations of p+ and p_ correspond respec- 
tively to first and second sounds. Now the right hand sides of Eqs. (17) and (18) 
should be VP and Vf respectively. We can write VP and Vf as functional of p+ 
and p_ by 

VP = AVp+-|-BVp_, (19) 
Vf = CVp+-FDVp_, (20) 



where 

A Sn ,^ , ^ \ , S22 , ^ N , gl2 ^ 

c = ^(P.+P-)-§(P.-P-) + 5^P-. 



The wave equations derived from Eq.(15)-(18) are reduced to 



^p_=CV^p++DV^p_. 



p+=AVV++BVV_, 

2^ 1 nv72. 



(21) 
(22) 
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Considering the plane waves that p+ and p_ oscillate around the equilibrium val- 
ues and with the frequency ft) and the wave number k like 

p+ = p+ + 5p+exp[;'(fe • r - cot)], 
p_ = p° + 5p_exp[j(fc-r - ft)f)], 
the sound velocities are 

where c = ft)/|fc|. These are obtained from the dispersion relation of the Bogoli- 
ubov excitations^ '" in two-component BECs in the limit of long wavelength. 

In superfluid '^He first and second sounds are modes that p and a indepen- 
dently oscillate. However, Eq. (23) does not necessarily describe first and second 
sounds because two modes are mixed. We can find that p+ and p_ oscillate inde- 
pendently when B and C vanish in Eqs. (21) and (22). This conditions are reduced 
to 

P? = Pi 

gll _ 822 
2 2 ' 

Then sound velocities of the two modes are 

4=,2±^, (24) 
mi 1712 

where i = ^ g\\p'(/m\ = \Jg72P^Jm\ and p'^ = p^ = p". The mode of c\ is 

oscillation of p+, first sound, and the mode of is oscillation of p_, second 
sound. First sound velocity increases with ^12 and second sound velocity decreases 
with g\2 (Fig.l). When |^i2| > 8^ \/8i\822, c+ or c_ becomes imaginary so that 
the dynamical instability leads to the collapse or the phase separation in the two- 
component BECs'". 

4 Summary 

We formulated the two-fluid model for two-component BECs, starting from the 
coupled GP equations. This model well corresponds to the two-fluid model in su- 
perfluid ^He expect for the mutual friction term. We obtained the condition that 
two sound modes are independent of each other like first and second sounds in 
superfluid *He. Second sound has an important role to investigate quantum turbu- 
lence in superfluid ''He. We are interested in how second sound interacts with a 
vortex in two-component BECs. In the future, we should investigate "mutual fric- 
tion" induced by the interaction between vortices and the Bogoliubov excitations 
in two-component BECs. The details will be reported soon elsewhere. 
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Fig. 1 Velocity of first and second sounds as a function of g[2. The solid and dashed line refers 
to first sound c\ and second sound respectively. 



Acknowledgements M.T. acknowledges the support of a Grant-in-Aid for Scientific Research 
from JSPS (Grant No. 21340104). 



References 

1. P. L. Kapitza, Nature 141, 74 (1937). 

2. L. Tisza, Nature 141, 913 (1938). 

3. L. Landau, J. Phys. U.S.S.R. 5, 71 (1941). 

4. C. J. Gorter and J. H. Melink, Physica. 15, 285 (1949). 

5. I. M. Khalatnikov, An Introduction to the Theory of Superfluidity (Perseus 
Publishing, 1965). 

6. K. Kasamatsu, M. Tsubota, and M. Ueda, Int.J. Mod. Phys. B 19, 1835 (2005). 

7. H. Takeuchi, N. Suzuki, K. Kasamatsu, H. Saito, and M. Tsubota, Phys. Rev. 
B 81, 094517 (2010). 

8. K. Sasaki, N. Suzuki, D. Akamatsu, and H. Saito, Phys. Rev. A 80, 063611 
(2009). 

9. Yu. G. Gladush, A. M. Kamchatnov, Z. Shi, P G. Kevrekidis, D. J. 
Frantzeskakis, and B. A. Malomed, Phys. Rev. A 79, 033623 (2009). 

10. C. J. Pethick and H. Smith, BOSE-EINSTEIN CONDENSATION IN DI- 
LUTE GASES Second Edition (Cambridge University Press, Cambridge, 
2008). 



